Search for: All records

The ultimate regularity of quantum mechanics creates a tension with the assumption of classical chaos used in many of our pictures of chemical reaction dynamics. Out-of-time-order correlators (OTOCs) provide a quantum analog to the Lyapunov exponents that characterize classical chaotic motion. Maldacena, Shenker, and Stanford have suggested a fundamental quantum bound for the rate of information scrambling, which resembles a limit suggested by Herzfeld for chemical reaction rates. Here, we use OTOCs to study model reactions based on a double-well reaction coordinate coupled to anharmonic oscillators or to a continuum oscillator bath. Upon cooling, as one enters the tunneling regime where the reaction rate does not strongly depend on temperature, the quantum Lyapunov exponent can approach the scrambling bound and the effective reaction rate obtained from a population correlation function can approach the Herzfeld limit on reaction rates: Tunneling increases scrambling by expanding the state space available to the system. The coupling of a dissipative continuum bath to the reaction coordinate reduces the scrambling rate obtained from the early-time OTOC, thus making the scrambling bound harder to reach, in the same way that friction is known to lower the temperature at which thermally activated barrier crossing goes over to the low-temperature activationless tunneling regime. Thus, chemical reactions entering the tunneling regime can be information scramblers as powerful as the black holes to which the quantum Lyapunov exponent bound has usually been applied. 

This paper reports the release of PathSum, a new software suite of state-of-the-art path integral methods for studying the dynamics of single or extended systems coupled to harmonic environments. The package includes two modules, suitable for system–bath problems and extended systems comprising many coupled system–bath units, and is offered in C++ and Fortran implementations. The system–bath module offers the recently developed small matrix path integral (SMatPI) and the well-established iterative quasi-adiabatic propagator path integral (i-QuAPI) method for iteration of the reduced density matrix of the system. In the SMatPI module, the dynamics within the entanglement interval can be computed using QuAPI, the blip sum, time evolving matrix product operators, or the quantum–classical path integral method. These methods have distinct convergence characteristics and their combination allows a user to access a variety of regimes. The extended system module provides the user with two algorithms of the modular path integral method, applicable to quantum spin chains or excitonic molecular aggregates. An overview of the methods and code structure is provided, along with guidance on method selection and representative examples. 

Path integral simulations reveal the mechanistic pathway of inter-ring excitation energy transfer in photosynthetic bacteria. 

Excitation energy transfer (EET) is fundamental to many processes in chemical and biological systems and carries significant implications for the design of materials suitable for efficient solar energy harvest and transport. This review discusses the role of intramolecular vibrations on the dynamics of EET in nonbonded molecular aggregates of bacteriochlorophyll, a perylene bisimide, and a model system, based on insights obtained from fully quantum mechanical real-time path integral results for a Frenkel exciton Hamiltonian that includes all vibrational modes of each molecular unit at finite temperature. Generic trends, as well as features specific to the vibrational characteristics of the molecules, are identified. Weak exciton-vibration (EV) interaction leads to compact, near-Gaussian densities on each electronic state, whose peak follows primarily a classical trajectory on a torus, while noncompact densities and nonlinear peak evolution are observed with strong EV coupling. Interaction with many intramolecular modes and increasing aggregate size smear, shift, and damp these dynamical features. 

The process of excitation energy transfer (EET) in molecular aggregates is etched with the signatures of a multitude of electronic and vibrational time scales that often are extremely difficult to resolve. The effect of the motion associated with one molecular vibration on that of another is fundamental to the dynamics of EET. In this paper we present simple theoretical ideas along with fully quantum mechanical calculations to develop a comprehensive mechanistic picture of EET in terms of the time evolution of electronic-vibrational densities (EVD) in a perylene bisimide (PBI) dimer, where 28 intramolecular normal modes couple to the ground and excited electronic states of each molecule. The EVD motion exhibits a plethora of dynamical features, which impart physical justification for the composite effects observed in the EET dynamics. Weakly coupled vibrations lead to classical-like motion of the EVD center on each electronic state, while highly nontrivial EVD characteristics develop under moderate or strong exciton–vibration interaction, leading to the formation of split or crescent-shaped densities, as well as density retention that slows down energy transfer and creates new peaks in the electronic populations. Pronounced correlation effects are observed in two-mode projections of the EVD, as a consequence of indirect vibrational coupling between uncoupled normal modes induced by the electronic coupling. Such indirect coupling depends on the strength of exciton–vibration interactions as well as the frequency mismatch between the two modes and leaves nontrivial signatures in the electronic population dynamics. The collective effects of many vibrational modes cause a partial smearing of these features through dephasing.